Summary

Thanks to feedback from colleagues, we have identified some errors in the SEDAC Global Human Influence Index dataset. Notably, there are some island areas that are heavily populated or otherwise obviously impacted by humans that return a false HII value of zero, indicating no human impacts. Here, we have removed these values from our dataset and reanalyzed the data. The results are almost identical to those in the main text of the paper, with no major interpretations changing. That said, we feel it important to point out the deficiency in this underlying dataset for future users.

We advise that future work utilize other global human impact databases that are more appropriate for coral reefs, such as Andrello et al “Global Human Pressure on Reefs” metric (referred to here and in our manuscript as WCS score) and Cinner et al’s Human Gravity.

Details on HII and WCS Score are provided in the following tabs.

HII details and error identification

The Global Human Influence Index (HII) is a globally gridded dataset with 1 km resolution grid cells. It was generated in 2005 by the NASA Socioeconomic Data and Applications Center (SEDAC) through a collaboration between the Wildlife Conservation Society (WCS) and the Columbia University Center for International Earth Science Information Network.

HII website and data download

Citation: Wildlife Conservation Society - WCS, and Center for International Earth Science Information Network - CIESIN - Columbia University. 2005. Last of the Wild Project, Version 2, 2005 (LWP-2): Global Human Influence Index (HII) Dataset (Geographic). Palisades, NY: NASA Socioeconomic Data and Applications Center (SEDAC). https://doi.org/10.7927/H4BP00QC.

HII Metric Details

HII is a composite metric of human influence that includes the following layers:
-Human population pressure (via human population density/ sq. km). Values from 0-10 assigned based on population density per square km
- Human land use and infrastructure (via urban areas, nighttime lights and land use/land cover)
- Urban area: if inside an urban polygon = 10, if outside = 0
- Light: ranked scores from 0-10 based on nighttime light values
- Land use: Urban areas = 10, irrigated ag = 8, rain-fed ag = 3, other cover = 0
- Human access (via coastlines, roads, railroads, and navigable rivers)
- Influence of coastlines was given 0 influence if > 15km from a coastline and 4 influence if within 15 km of a coastline
- Roads were assigned an influence score based on proximity. If a site was within 2 km of a major road, it received the highest score (8), between 2 and 15 km received a 4, and > 15 km from a major road received a 0
- Railroads were scored as 0 if > 2 km from a railroad and 8 if within 2 km of a railroad - Navigable rivers received a 0 score if > 15 km from a river and 4 if within 15 km of a river
The metric is calculated by assigning an influence score to various values in each layer (higher = more human influence). After values are assigned to each layer of the metric, they are summed to give an HII value.

HII metric calculation website

Error Identification

Through expert feedback and detailed exploration of the dataset, we have identified a number of location in which HII is falsely reported to by zero (HII=0). An HII value of zero suggests no human impact on that site.

Figure 1: Map of HII data from the SEDAC webpage. Gray areas indicate missing data. In these locations, HII=0.

As a result of this missing data in the underlying HII layer, we combed our dataset and assessed each site in which HII=0 using Google Earth. Since we calculated HII for the manuscript as the sum of all HII within 100km of each site, we created a 100km radius circle on Google Earth with the lat/long of the site at the center and assessed (visually) whether there was clearly human influence or not. Since HII was generated in 2009, we utilized modern satellite imagery (the most recent available images) and images from 2008 and/or 2009 at each site. In some cases, there was very clearly human influence (Fig 2). In these instance, we removed the point from our dataset and reanalyzed. In other cases, there was no visible human influence within 100km. These sites were retained (Fig 3).

Fig 2: Maduvvari had clear human influence and was removed from the dataset.

Figure 3: Cocos Island, Costa Rica had no visible human influence and was retained.

Data Removal

After exploring the data, we removed 41 sites from our Recovery database (out of an original 182) and 40 sites from our resistance database (out of an original 184). Most of these removals were in the Indian Ocean. It seems that the missing data in the HII dataset are concentrated in small island states, specifically within the Indian Ocean. We urge caution in using this dataset for this region of the world in the future and note that such gaps do not appear to exist in the Andrello et al WCS score or the Cinner et al Human Gravity data.

WCS Score Details

Explore these tabs for details on the WCS Score metric (along with links to download the data)

WCS Score Background

WCS Score is derived from a pre-print written by Marco Andrello, Emily Darling, Amelia Wenger, Andres F. Suarez-Casto, Sharla Gelfand, and Gabby Ahmadia that is available on BiorXiv. The paper is a collaboration between The National Research Council (Italy), MARBEC (University of Montpelier, France), Wildlife Conservation Society (WCS), and a number of other institutions. Link to paper

The Supplementary Info are very helpful and contain details of each component of the metric Download Supplementary Info

Data and Code

Citation: A global map of human pressures on tropical coral reefs Marco Andrello, Emily Darling, Amelia Wenger, Andrés F. Suárez-Castro, Sharla Gelfand, Gabby N. Ahmadia bioRxiv 2021.04.03.438313; doi: https://doi.org/10.1101/2021.04.03.438313

Metric Details

‘WCS Score’ is a composite metric that utilizes the following layers:
- Tropical coral reef locations (from Beyer et al 2018, based on UNEP-WCMC dataset from 2010). Resolution is ~ 5km.
- Fishing (Market Gravity from Cinner et al 2016, 2018, resolution = 10km)
- Coastal Development (gridded population of the world from SEDAC, 2018 – an updated version of the same data used in HII. They used number of people living within a 5km buffer of each reef polygon as the actual metric here)
- Industrial Development (World port locations- ports were identified from a google data product, given a point location, and then all port points within 5 square km of a reef were summed)
- Tourism (This value is directly from Spalding et al 2017 – an estimate of the annual economic value of each reef)
- Sediment Pollution (Via a newly calculate metric for this paper- Sediment delivery from rivers was estimated in 2 ways based on land cover, management, rainfall, discharge, land use, and more. Sediment exposure to reefs was modelled using sediment plume models and estimated at tons of sediment per square km)
- Nitrogen Pollution (Using FAO data on catchment level crop cover and national nitrogen fertilizer use, they modeled nitrogen plumes from rivers to estimate nitrogen pollution)

Each metric was assigned a weight (1-6 with 6 meaning most impactful) for each reef pixel using perceptions of risk from a study by Stephanie Wear (2016). A weighed average of the 6 metrics was used to generate a cumulative score, or as we are calling it, ‘WCS Score’ for each pixel.

See supplementary methods from Andrello et al for more detail on all of these metrics and how they were calculated.

Updated Tables and Figures

Tables

Updated Recovery Model Table
Parameter Coefficient SE CI CI_low CI_high t df_error p Effects Group
(Intercept) 0.6159553 1.4456038 0.95 -2.2520815 3.4839921 0.4260886 100 0.6709590 fixed
regionE. Pacific 0.8626023 1.6261043 0.95 -2.3635423 4.0887468 0.5304717 100 0.5969610 fixed
regionIndian Ocean 0.4907281 1.6541673 0.95 -2.7910928 3.7725490 0.2966617 100 0.7673401 fixed
regionW. Pacific 2.2685120 1.4381334 0.95 -0.5847037 5.1217276 1.5774003 100 0.1178618 fixed
disturbanceBleaching, COTS -1.0613801 2.0570847 0.95 -5.1425776 3.0198174 -0.5159633 100 0.6070190 fixed
disturbanceBleaching, Disease -7.1914099 2.9729602 0.95 -13.0896783 -1.2931415 -2.4189392 100 0.0173734 fixed
disturbanceBleaching, Storm -3.3884743 1.4679047 0.95 -6.3007553 -0.4761932 -2.3083749 100 0.0230372 fixed
disturbanceCOTS -1.5277335 1.2577732 0.95 -4.0231198 0.9676528 -1.2146335 100 0.2273663 fixed
disturbanceDisease -2.4312290 2.0111592 0.95 -6.4213116 1.5588536 -1.2088695 100 0.2295636 fixed
disturbanceStorm 0.8273779 1.2119434 0.95 -1.5770832 3.2318391 0.6826870 100 0.4963831 fixed
hii100km2 0.5509937 0.4735701 0.95 -0.3885558 1.4905433 1.1634894 100 0.2473997 fixed
MPA_during_recovYes -0.1571351 0.6814394 0.95 -1.5090915 1.1948212 -0.2305930 100 0.8181021 fixed
recovery_time2 -0.4946426 0.4063243 0.95 -1.3007784 0.3114932 -1.2173592 100 0.2263326 fixed
distance_to_shore_m2 -0.0644055 0.4004037 0.95 -0.8587951 0.7299841 -0.1608515 100 0.8725349 fixed
cml_scr2 0.1607983 0.4121551 0.95 -0.6569057 0.9785023 0.3901402 100 0.6972626 fixed
pre.dist.cc 0.0252404 0.0174126 0.95 -0.0093057 0.0597865 1.4495491 100 0.1503132 fixed
SD (Intercept) 2.4420087 NA 0.95 NA NA NA NA NA random study
SD (Observations) 2.0059320 NA 0.95 NA NA NA NA NA random Residual
Updated Resistance Model Table
Parameter Coefficient SE CI CI_low CI_high t df_error p Effects Group
(Intercept) 11.5106568 4.7955428 0.95 1.9690309 21.0522827 2.4002824 81 0.0186757 fixed
regionE. Pacific -5.7378136 7.2544007 0.95 -20.1717955 8.6961683 -0.7909425 81 0.4312874 fixed
regionIndian Ocean 1.7274082 6.1714211 0.95 -10.5517840 14.0066004 0.2799044 81 0.7802644 fixed
regionW. Pacific -7.7074993 5.1701283 0.95 -17.9944328 2.5794342 -1.4907753 81 0.1399046 fixed
disturbanceBleaching, COTS -10.5828188 7.3641010 0.95 -25.2350699 4.0694323 -1.4370822 81 0.1545463 fixed
disturbanceBleaching, Storm -12.4616494 6.6555358 0.95 -25.7040780 0.7807792 -1.8723736 81 0.0647635 fixed
disturbanceCOTS -10.4173519 4.7811930 0.95 -19.9304263 -0.9042775 -2.1788185 81 0.0322513 fixed
disturbanceDisease -8.5517700 8.7566475 0.95 -25.9747517 8.8712117 -0.9766032 81 0.3316734 fixed
disturbanceStorm -6.5298251 4.3023782 0.95 -15.0902082 2.0305580 -1.5177246 81 0.1329773 fixed
hii100km2 -1.8204206 2.2172273 0.95 -6.2320074 2.5911662 -0.8210347 81 0.4140364 fixed
MPA_during_resistYes -5.4074026 2.9658521 0.95 -11.3085179 0.4937127 -1.8232206 81 0.0719592 fixed
resistance_time_2 -2.4712555 1.8113986 0.95 -6.0753705 1.1328596 -1.3642803 81 0.1762580 fixed
cml_scr2 -1.5865288 1.8485066 0.95 -5.2644770 2.0914194 -0.8582760 81 0.3932727 fixed
travel_time.tt_pop2 -5.0241856 2.0619613 0.95 -9.1268417 -0.9215294 -2.4366052 81 0.0170182 fixed
pre.dist.cc -0.6933835 0.0694121 0.95 -0.8314918 -0.5552753 -9.9893781 81 0.0000000 fixed
SD (Intercept) 9.1492581 NA 0.95 NA NA NA NA NA random study
SD (Observations) 7.0576701 NA 0.95 NA NA NA NA NA random Residual

Figures

Updated Figure 1 (regional trends) Figure 1 is very similar to the published figure. The only noticeable difference is in the Indian Ocean histrogram in panel A (Recovery Rate), where the distribution is unimodal instead of bimodal.

Updated Manuscript Figure 2 (model coefficient plots) Figure 2 has changed slightly. For this version, WCS SCORE is called ‘cml_scr2’. There are slight changes to some of the model coefficients following the removal of false HII=0 sites. Namely, HII is no longer (significantly) positively correlated with Recovery Rate, though that trend still appears to be present in Figure 3.

Updated Manunscript Figure 3 (HII/WCS Score vs. Recovery/Resistance) Figure 3 has changed slightly from the main text. Panels A and B there are fewer HII=0 points, as many of these were false zeros and have been removed in the re-analysis. There were no apparent trends in Panels B, C, and D in the orignal text and this remains true. In Panel A, though there is a hint of a positive correlation between HII and Recovery Rate, this is no longer statistically significant. In short, there are no clear relationships between either human influence metric and Recovery or Resistance.